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We systematically improve the recent variational calculation of the 
imaginary part of the ground state energy of the quartic anharmonic 
oscillator. The results are extremely accurate as demonstrated by deri- 
ving, from the calculated imaginary part, all perturbation coefficients 
via a dispersion relation and reproducing the exact values with a re- 
lative error of less than fO~^. A comparison is also made with results 
of a Schrodinger calculation based on the complex rotation method. 

1) It is well known that the evaluation of perturbation expansions can be 
greatly improved with the help of a variational approach [1]. Recently it 
was found that also the calculation of tunneling processes can be carried 
far beyond previously available semiclassical approximation [2] by combining 
it with variational methods [3]. The range of validity of the semiclassical 
results was extended far into the regime of low barriers and large tunneling 
rates. Variational methods allow us to calculate decay rates even when the 
tunneling barriers are so low that they do not contain any bound states at 
all. These are the sliding rates of Ref. [3]. For the quantum mechanical 
anharmonic oscillator with a potential 



low-order calculations of tunneling and sliding rates yield imaginary parts 
which are so accurate that a dispersion relation produces perturbation coef- 
ficients which agree with the exact values calculated by Bender and Wu [4] 
to all orders in perturbation theory. 

In this note we want to investigate the increasing accuracy of the ap- 
proximation method when carrying the systematic improvement scheme of 
ref. [1] to increasingly higher orders. At the end we shall obtain the imaginary 
part of the ground state energy with a precision of better than 0.01 %. 
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Up to now, the best developed theory of tunnehng processes apphes only 
to the limit of high potential barriers or small values of —g where a semiclas- 
sical treatment is appropriate. For increasing — </, the fluctuation corrections 
to the semiclassical results become increasingly difficult to calculate. For 
large — </, where the decay proceeds by sliding rather than tunneling, no field 
theoretic method was available until [3]. The power of the variational me- 
thod is to obtain accurate results over the entire negative (/-axis, including 
the two difficult regimes. 

2) The sliding regime was found to be easily calculable. It merely requires 
an analytic continuation of the variational approximations which is known 
for positive g. There the Rayleigh-Schrodinger perturbation theory yields a 
power series expansion ifi = 1) 

where ef ™ are rational numbers 
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The superscript BW refers to Bender and Wu [4] who have found a recursion- 
relation by which one can easily find the first few hundred terms. The dis- 
advantage of the series (2) is that it has a zero radius of convergence due to 
the factorial growth of the coefficients ef™. To be evaluated approximately 
it has to be truncated at a finite order. The best results are obtained by 
going to an order 3/4(/. Reliable numerical values cannot be expected for 
g>0.1. 

The evaluation is greatly improved by a variational treatment. It proceeds 
as follows: The harmonic term of the potential is split into an arbitrary 
harmonic term and a remainder: 
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V{x) = ^x' + (5) 
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with an interaction 

4 g 
one performs perturbation expansion in powers of at a fixed r: 

i?,(5,r) = 0^e,(r)(^) (7) 

The calculation of the new series up to a specific order k does not require 
much additional work since it is easily obtained from the ordinary perturba- 
tion series (2) by replacing u by + gr l2 and expanding in powers of g 
up to the A;th order. This yields 

The truncated power series 

Wk{g,^):=Ek{g,2{u^ -9?)lg) (9) 

is certainly independent of in the limit of large k. At any finite order k, 
however, it does depend on 0, the approximation having its fastest speed 
of convergence where it depends least on 0. Denote the order-dependent 
optimal value of by ilk{g). The quantity Wk{g^ ftk{g)) is the new approxi- 
mation to E{g) [5]. Since Wk{g,ft) is of the form Em=o with 
polynomials am{g)^ the defining equation for its extrema is a polynomial of 
order 3k in 0. Thus only in the case of odd A;, a minimum is guaranteed to 
exist. In the case of even A;, the alternating sign of ao[g) turns out to remove 
the minimum for smaller k. Not considering very small positive g (for which 
there always exist k real solutions in the vicinity of io) it has been found 
numerically that VK/;((/,0) has no extrema for k = 2,4, one extremum for 
A; = 1, 3, 5, 7, 9, 11, 13, two extrema for A; = 6, 8, 10, 12, 14, 16, 18, 20 and three 
extrema for k = 15,17,19,21. At further increasing orders the situation 
becomes increasingly difficult to handle numerically because the numerous 
extrema are to close. An interesting situation occurs in the case of three 
extrema when there are two minima. The global minimum of Wk{g^ 0) turns 
out to yield a 100-times less accurate energy value than the higher lying local 
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minimum. In Table 1 we display the values obtained from the better local 
minimum in the cases k = 15 and k = 21. In this way, an agreement of 
0.01 ppm with the values given in the literature [6] was reached. An inte- 
resting fact is that the inequality -Eexactl*/) ^ Min 0)) which in the 
case k = 1 follows from the Jensen-Peierls inequality seems to hold also for 
A; = 3, whereas it is definitely violated for higher k. Note, however that the 
above-discussed preferable local minima do again provide upper bounds in 
accordance with the Jensen-Peierls inequality. 

Our approximation is easily generalized to the excited states. One only has 
to replace the ground state perturbation coefficients by the corresponding 
coefficients ef^ for the nth excited state. We have calculated them up to 
k = 5 for all n using usual Rayleigh-Schrodinger perturbation theory. The 
error of the resulting approximation behaves similar to that of the ground 
state energies. 
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0.5591463 
0.6379918 
0.6961758 
0.8037706 
0.9515685 
2.4997088 
3.9309313 
6.6942208 
13.366908 
18.137229 



0.5603074 
0.6416299 
0.7016616 
0.8125000 
0.9644036 
2.5475804 
4.0084609 
6.8279533 
13.635283 
18.501659 



0.5591542 
0.6380358 
0.6962536 
0.8039140 
0.9517998 
2.5006996 
3.9325538 
6.6970329 
13.372561 
18.144908 



0.5591462 
0.6379899 
0.6961717 
0.8037615 
0.9515517 
2.4996213 
3.9307858 
6.6939668 
13.366395 
18.136533 



0.5591462 
0.6379897 
0.6961712 
0.8037601 
0.9515490 
2.4996061 
3.9307603 
6.6939221 
13.366305 
18.136411 



0.5591462 
0.6379903 
0.6961723 
0.8037623 
0.9515526 
2.4996219 
3.9307863 
6.6939673 
13.366396 
18.136534 



0.5591463 
0.6379918 
0.6961758 
0.8037707 
0.9515685 
2.4997089 
3.9309316 
6.6942213 
13.366908 
18.137230 



0.5591463 
0.6379918 
0.6961758 
0.8037706 
0.9515685 
2.4997088 
3.9309313 
6.6942209 
13.366908 
18.137229 



Table 1: Comparison of the exact values of E[g) with those obtained by 
minimization of Wk{g^ 0). 



To judge the quality of our approximation consider the asymptotic 
strong-coupling behavior of E[g) both for negative as for positive g which's 
exact numerical values are known with great precision [6]. It can be para- 
metrized as follows: 

E.^{g) = {glAfI' [a^ + (i^{glA)-'" + inigl^Y'" • • •] • (10) 

The asymptotic behavior of Im£'„( — </) is obtained by an analytic continua- 
tion: 

l^E^{-g) = {gjif' [a; + fUdl^r"' + I'Aal^r''' • • •] , (11) 
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/3o 


k 


"0 


/3o 


1 


0.681420222 


0.13758 


9 


0.667958214 


0.143700 


3 


0.668269437 


0.14347 


11 


0.667960613 


0.143698 


5 


0.667960581 


0.14369 


13 


0.667962642 


0.143696 


7 


0.667956063 


0.14370 


15 


0.667986302 


0.143669 








exact 


0.667986259 


0.14367 



Table 2: Comparison of the exact asymptotic behavior with those obtained 
by our approximation at increasing order k of the variational approximation. 

with coefficients 

< = ^0!n , f^'n = -^f^n , 7n = 'In ■ (12) 

Table 2 shows the values of our approximation for the ground state coeffi- 
cients cto and /3o in comparison with the exact values. The agreement is seen 
to improve rapidly with the order k of the approximation. Our approxima- 
tion for the coefficients a„ and /3„ of the excited states is displayed in Table 
3. The order of the approximation is A; = 5. 

3) It was shown in ref. [3] for the lowest case A; = 1, that the above approxima- 
tion scheme makes also sense for sufficiently negative g. There the imaginary 
part of the ground state energy is obtained by an analytic continuation of the 
minimum of Wi{g,Vt) = 3^/160^ + 1/40 + 0/4. For g < -ij'i^l'^ ^ -0.25 
the optimal value of moves into the complex (/-plane and the imaginary 
part of the energy describes a sliding process into the abysses. In the case 
of higher (odd) orders A;, this analytic continuation has to be done numeri- 
cally. There are standard computer programs which give all the 3k roots of 
JVKfc((/, 0)/ JO for arbitrary complex g to any precision. From these we have 
to choose the one that is continuously connected (i.e. on a semicircle in the 
complex g-plane) with the best (in the sense as discussed above) solution 
on the positive (/-axis. With increasing orders the range of validity of the 
approximation moves closer and closer to the origin. This is quite amazing 
since the behavior of lniE[g) in the limit of small —g possesses an essential 
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1 

2 
3 
4 
5 
6 
7 
8 
9 

10 



0.6679605808 
2.3922963623 
4.6688449582 
7.2862672627 
10.171684571 
13.282179566 
16.588882965 
20.070767345 
23.711695355 
27.498813862 
31.421596541 











1 
1 
1 
1 



1436962401 
3586118166 
5059561947 
6340450238 
7501258899 
8577545339 
9589691819 
0550754253 
1469734787 
2353174155 
3206014789 



Table 3: Coefficients governing the asymptotic behavior of the excited states 
at a fixed order A; = 5 of the variational approximation. 

singularity 



Our analytic approximants approach rapidly this nonanalytic behavior. Fi- 
gure 1 shows the approximations at different orders. We have plotted the 
reduced imaginary part e{g) = lniE[g) /lniEsc{g) whose semiclassical limit is 
unity. 

To check the correctness of our result we compare it with the values given 
in the literature which have an accuracy of 0.1 % [7]. We derive even more 
precise values by using the so-called "complex rotation method" [10]. It yields 
a matrix recurrence relation which can be solved using matrix continued 
fractions (see ref. [11]). This method yields numerical values which coincide 
with those of our approximation within a relative error of 0.1%. Table 4 
compares the three independently obtained results. Our values have about 
the same precision as those obtained by Drummond. This can be seen from 
the results of the "complex rotation method" which are calculated to very 
high accuracy. 




(13) 
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Figure 1: Reduced imaginary part at different orders of our approximation 
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4332072064 





4333 


0.4333195993 


-1.4 





4534503231 
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0.4535431192 


-1.6 





4763745530 





4764 


0.4764405650 


-1.8 





5009071563 





5009 


0.5009448702 


-2.0 





5264438875 





5265 


0.5264541474 



Table 4: Comparison of our approximation of the reduced imaginary part 
e{g) (order k = 15) with Drummond's results and the precise numerical 
values obtained via the solution of the Schrodinger equation with the complex 
rotation method. 
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4) In the vicinity of the origin where the decay proceeds by tunnehng, the 
previous approximation gives no imaginary part. But we can adjust it to the 
semiclassical behavior as follows: One first has to calculate corrections to the 
semiclassical behavior 

lniE[g) = Im£'^((/) exp 

where dj are positive rational numbers. The first of them can be calcula- 
ted using field theoretic techniques [9]. In the special case of the quartic 
anharmonicity Zinn- Justin [8] has found them numerically up to diQ. Unfor- 
tunately, the power series in g is only asymptotic one so its accuracy does not 
increase if k is carried beyond 4/3|(/|. The variational approach improves 
the results and the convergence dramatically since the effective coupling con- 
stant is the dimensionless quantity a[g) = g/Q,^[g) which is always smaller 
than 0.6 this value being reached only in the strong-coupling limit — )• oo 
(this limit holds for k = 1; for k = 3 the maximal value decreases to 0.3. For 
a plot see [12]). 

To improve this variationally one rewrites this by introducing a trial 
frequency as 

lniE[g) = lniE2.[g) exp 

with 

J_i = 1, do = 0, ^1 = (16) 

At this point we replace co by _|_ g/2r , ( r = 2[(jj'^ — O^)/*/ ) and 

expand in powers of g at fixed r up to order k which yields 

lniE[g) = Im£'g^((/) exp 

where 

m=0 \ / 



k 



3|g| 
"4cu3 



(14) 



^iog(^)- y d, 
j=-i 



3|g| 
"4cu3 



(15) 



02 



02 



j=-i 



(17) 
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is still determined by an analytic continuation of the above optimal Oi for 
all k. For k = 1 we obtain therefore 



lmE{g) = lmE^^{g) exp 

where we have used the abbreviation x = (cu^ — 0^)/0^. The resulting ap- 
proximation is shown in Figure 2 for k = 1 and k = 2. 

The same reasoning can also be used for the excited states. Difficulties 
arise only in the calculation of dk^n for arbitrary n since the field theoretic 
method of [9] cannot be applied. There exist however a WKB calculation 
of di^n [4]. Figure 3 shows these variational approximations for the ground 
state and the first two excited states. We have also plotted the usual appro- 
ximation for k = 1 and A; = 3 to show how well the two approximations fit 
together. Figure 1 and 2 are magnifications of the dashed box. 

We can now use our knowledge of the entire cut to calculate the pertur- 
bation coefficients ef™ by inserting it into the dispersion relation 

er = - r lmE{g)g-'-Ug. (20) 

TT J-oo 

If we take for example k = 15 we can use our approximation left of —2/15 ~ 
0.133 and Fq. (14) truncated after A; = 10 to continue it up to = 0. The 
results are given in Table 5 together with the values for k = 3 (which were 
calculated in combination with the semiclassical result as in ref. [3]). They 
show again the quality of our approximation. 

5) The calculations presented in this paper were only intended as a prepara- 
tion for finding the imaginary part of the correlation functions of quantum 
field theories. If we succeed to extend our method to the (^"^-theory we would 
be able to study its properties with a much higher accuracy than presently 
available on the basis of standard resummation methods. In particular, there 
is hope that the precision of renormalization group calculations of critical in- 
dices can greatly be improved. 
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(19) 
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Figure 2: Approximation for the reduced imaginary part in the tunnehng 
regime 
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Figure 3: Approximation of the reduced imaginary part of the ground state 
and the first two excited states 
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k = 15 
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0.75 


0.749321679 


0.750035087 
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2.625 


2.614620115 


2.624986293 
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20.8125 


20.7186128 


20.81251773 
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241.289063 


240.857317 


241.289495 
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3580.98047 


3590.69587 


3580.97722 
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63982.8135 


64432.53875 


63982.3799 
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1329733.73 


1342857.03 


1329709.29 


8 


31448214.7 


31791078.07 


31447170.8 


9 


833541603 


842273537 


833504391 


10 


24478940700 


24703889150 


24477781497 



Table 5: The exact perturbation coefficients and those calculated from our 
approximation and the dispersion relation. 
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